Noise Correlations in one-dimensional systems of ultra-cold fermions 
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Time of flight images reflect the momentum distribution of the atoms in the trap, but the spatial 
noise in the image holds information on more subtle correlations. Using Bosonization, we study such 
noise correlations in generic one dimensional systems of ultra cold fermions. Specifically, we show 
how pairing as well as spin and charge density wave correlations may be identified and extracted from 
the time of fiight images. These incipient orders manifest themselves as power law singularities in the 
noise correlations, that depend on the Luttinger parameters, which suggests a general experimental 
technique to obtain them. 



Recent advances in trapping and control of ultra cold 
atoms prompted the realization of truly one dimensional 
systems. Important benchmarks, such as the Tonks- 
Girardeau eas [J, 0] and the Mott transition in one 
dimension|3|, have been achieved by trapping bosonic 
atoms in tight tubes formed by an optical lattice poten- 
tial. Novel transport properties of one dimensional lattice 
bosons have been studied using these techniques [4] . More 
recently, a strongly interacting one dimensional Fermi 
gas was realized using similar trapping methods ^5]. In- 
teractions between the fermion atoms were controlled by 
tuning a Feshbach resonance in these experiments. 

Even the simplest model of interacting spinless 
fermions in one dimension exhibits different regimes of 
behavior. Although spinless fermions confined to an iso- 
lated line display a single thermodynamic phase, it is 
convenient to characterise different parameter regimes in 
terms of which susceptibilities are divergent. A divergent 
susceptibility in a particular channel implies that cou- 
pling independent one dimensional systems in that chan- 
nel will lead to an instability. If we classify 'phases' of the 
one dimensional system in this way, then the 'phase' dia- 
gram of spinless fermions includes a diverging charge den- 
sity wave (CDW) susceptibility for repulsive interaction 
and diverging superconducting (SC) susceptibility for at- 
tractive interaction. In the case of spin-1/2 Fermions, 
which can exist either in a spin gapped or gapless phase 
in the continuum, a much richer competition of 'phases' 
(as defined from divergent susceptibilities) occurs. These 
include both triplet SC (TSC) and singlet SC (SSC)as 
well as CDW and spin density wave (SDW) (see e.g. [a]). 
Numerous proposals were given for realization of a vari- 
ety of different phases in ultra cold Fermi systems 
as well as Bose-Fermi mixtures |3, fojl. However methods 
for detection and experimental characterization of all of 
these strongly correlated states are lacking. 

There is a fundamental difficulty, not confined to 
atomic systems, in applying a unified probe of correla- 
tions. A standard detection scheme in solid state systems 
involves coherent scattering of external particles, such as 
neutrons or photons, that naturally couple to the sys- 



tem (weakly) via the particle hole channel. Such a probe 
is therefore sensitive to spin and charge order, but not 
to superconductivity. Thus parts of the phase diagram 
of one dimensional systems, which in theory are dual to 
each other, are not treated equally by such probes. In 
this paper we show that by measuring the noise correla- 
tions in time of flight images one can gain access to, and 
distinguish between different power-law particle-particle 
correlations, such as TSC and SSC, and particlc-hole cor- 
relations such as CDW and SDW. 

Atomic shot noise in time of flight imaging was pro- 
posed in Ref. [lo| as a probe of many body correlations 
in systems of ultra cold atoms. Speciflcally, what may 
be measured using this approach is the momentum space 
correlation function of the atoms in the trap 
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where a, a' are spin indices. The method was success- 
fully applied to experiments measuring correlations in 
both particle-particle and particle hole channels. For ex- 
ample Greiner et al[ll] measured pairing correlations in 
an expanding cloud of atoms originating from dissociated 
molecules. Other experiments measured density correla- 
tions in insulating states of lattice bosons [12, ll3| and 
fermions [l^. A related effect corresponding to the clas- 
sical limit of the noise correlations was observed in Ref. 

The analysis In the examples considered in [10(, cal- 
culation of the noise correlations ^ was greatly simpli- 
fled by the assumption of a mean held BCS state and 
a perfect Mott insulator respectively. However, mean 
fleld theory fails to describe the one dimensional systems 
discussed here. Instead the long distance behavior of cor- 
relations is described by the Luttinger liquid framework, 
which yields power law decays, determined by one or two 
Luttinger parameters. We shall calculate the correlation 
function ^ within this framework, and extract its lead- 
ing singularities. These turn out to be directly related to 
the pure CDW, SDW and SC correlations. 

Spinless fermions - Let us start by treating in detail 
the simplest case, namely that of spinless fermions in the 
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FIG. 1: The noise correlation function (|6]) of a spinless Lut- 
tinger liquid for values of the Luttinger parameter K a,) - 
d): K = 0.4, 0.8, 1.25 and 2.5. For K < 1/2 there are strong 
anti correlations (reflecting the particle-hole nature of CDW) , 
which appear as a power law singularity on the line q — q' — 0. 
For K > 2 correlations switch to particle-particle (SC) and 
appear as a powerlaw singularity on the line q + q' = 0. For 
1/2 < K < 2 we flnd a singularity only when the two mo- 
menta approach the two fermi points (i.e g = g' = 0). The 
system shows a signature of both particle-hole and particle- 
particle correlations with the sign of correlation depending on 
the directions of approach to the fermi points. 



continuum. Later we shall generalize to spin-i fermions. 
To implement a Luttinger liquid of spinless fermions one 
can use fully polarized fermionic atoms. The interaction 
may be tuned via Feshbach resonance in an odd angular 
momentum channel (e.g. p-wave). Another possibility 
for tuning interactions is to use a Bose-Fermi mixture, 
where phonons in the bose superfluid mediate effective 
interactions between the fermions 0. 

Regardless of the microscopic interactions, the low en- 
ergy physics of one dimensional spinless fermions is de- 
scribed by a universal harmonic theory, the Luttinger 
liquid. The action, describing phonons in the Fermi sys- 
tem, is written in terms of the phase field $ or its dual 
8: 
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Here fj, — 0,1, xq = vpr and K is the Luttinger parame- 
ter. K > 1 {K < 1) corresponds to attractive (repulsive) 
interactions. To calculate physical correlation functions 
one needs the bosonization identity connecting the oscil- 
lator fields to the physical fermions [lit • 
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where a is an infinitesimal to be taken to zero at the 
end of any calculation. 'ijj±{x) correspond to the slowly 



varying, left and right moving components of the fermion 
field ij{x) = e~'-''f'^ip-{x) + e'''"f'^'ip+{x). The commuta- 
tion relation of the oscillator fields is given by: 
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The relevant order parameters can be expressed sim- 
ply in terms of the slow fields. The 2k f charge 
density wave operator order parameter is given by 
OcDwi^) — e'^^'^^^ipl.ip-, while the SC order parameter 
is Oscix) = ip+Tp-. A standard calculation fij using ^ 
and ^ yields {Ol;jjy^{x)0 CDw{^)) ^ cos(2fcj^)|a;| 
and (Ot;c(a;)Osc(0)) ^ \x\-'^/'^. Taking the spatial and 
temporal Fourier transform one can see that CDW and 
SC susceptibilities diverge when K < \ and K > 1 re- 
spectively. 

We are interested in the equal time correlation function 
([1]) with k and k' close to the opposite fermi points at 
±fc/. The operator {rikUk'), with k (k') near the fermi 
points kf {—kf) is given by: 
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where q — k—kf and q' = k' +kf denote the momenta rel- 
ative to the respective fermi points. Using the bosoniza- 
tion identity ^ and the commutation relation ^ with 
the harmonic action ([2]) we obtain for the noise correla- 
tion function ([1]): 

/piqxi2+iq' X34 
{2^YV •^(^i2)-^*(^34)(^-l).(6) 
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the integration is over X12, 0:34, 0:23, 
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with a a short distance cutoff (which will be kept finite 
in contrast to a, see e.g. @), and exponents g = {K + 
K-^ - 2)/4 and h = {K - K-^)/A. We shaU undertake 
a full numerical integration of ([6]) . But first let us obtain 
the leading singularities in the limits if <C 1, K ^ 1, and 
|1 - i^l < 1. We note that JS]) 'contains' both the CDW 
and SC correlations. For example, the contribution to 
([6]) from setting xi ^ x^ = x and X2 ~ x^ = x' is simply 
the fourier transform of the CDW correlation at q — q', 
with a minus sign due to anticommutation of the Fermi 
fields. For K ^ 1 since CDW correlations decay very 
slowly, we expect that there are no singular contributions 
arising from the deviations xn, X23, so that the CDW 
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correlations dominate the integral. Indeed, xn and 0:23 
appear in the numerator of A, and since for if <C 1, 
h is very negative, the contribution of large deviations 
is effectively suppressed. This gives G{q,q') cc —\q — 
qi]2K-i ^ In a similar manner setting xi — xj, and X2 = x^ 
yields the SC correlation at q + q' . For K ^ 1, h is 
large and positive, effectively suppressing the deviations 
Xi3 and 3:24. Thus in this limit we find G{q,q') oc |g + 
g'|2/-ff-i With these considerations we would not expect 
a singularity in G{k, k') in the region ^ < K < 2 due to 
SC or CDW correlations. However, the approximation 
presented above is not valid for such moderate values of 
K. Within this regime we can obtain analytic results for 
the leading singularities of ([6]) in the limit — 1| <C 1 
{\h\ <^ 1). The idea is to expand G{q, q') to powers of h. 
The term linear in h vanishes identically, and to quadratic 
order: 



G{q, q') ~ — /i^sgn((7)sgn(g')min 
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We now turn to integrate ^ numerically in order to 
obtain G{q,q'). This is carried out on a finite size sys- 
tem with periodic boundary conditions so that xtj in ^ 
are replaced by (L/27r) sin(27rxij/L). The results, sum- 
marized in Fig. [1] indicate that the anticipated power 
law singularities persist in the entire region K > 2 and 
K < 1/2, where the equal time order parameter correla- 
tions yield a singular contribution. The structure of the 
correlations in the regime 1/2 < X < 2 is consistent with 
the analysis presented above for K 1. 

To predict the visibility of the peaks we consider the 
normalized correlation function typically used in such 
experiments 
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The numerator is given by the connected part of the 
correlation function ([5]), i.e. the part containing A, 
while the denominator is the disconnected part. Fig. [2] 
shows a computation of G{k) for a realistic system of 
decoupled tubes with 100 particles per tube. Thus the 
short range cut-off, which appears in Eqs. ([7]) is set to 
a = 2'K/kf = L/100. The momenta are then given in 
units of 27r/L, which appears as the scale on the hori- 
zontal axis. The presence of QLRO, rather than LRO, 
is evident by the appearance of broad power-law dip in 
C(k) compared to the sharp dip in the case of perfect 
order. We note that the distinct dip is clearly observable 
only for strong interactions (i.e. K < 0.25). 

Spin-^ fermions - A spin-^ system exhibits much 
richer physics than the spinless one, and it may be re- 
alized using interactions in the s-wave channel. The low 
energy physics is captured by the universal action: 
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FIG. 2: The normalized correlation function (|9j for a real- 
istic finite system with N « 100 particles per Id tube, as a 
function of momentum k in the vicinity of 2k f. (The scale 
of the horizontal axis is indicated by the marks which cor- 
respond to increments of 2tt/L.) The result is plotted for 
K = 0.3,0.2,0.18 in the CDW regime. The broad dip due 
to the power-law singularity should be contrasted with the 
sharper dip for the case of a finite Mott insulating state (MI) . 
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Where v ~ p,a correspond to the charge and spin de- 
grees of freedom, and Qp^a = (©T ©i)/v^ and ^p,^ = 
($1 ± $1)/ \/2. The nonlinear term that arises in the spin 
channel is due to back scattering. The bosonization for- 
mula is identical to ([3]) with Fermi as well as oscillator 
fields taking a spin index. 

The low energy fluctuations of the spin-i LL are domi- 
nated by four order parameters constructed from the slow 
fermion fields: Ocdw = Yjs^'^'^'''i'-,s'4'+.s, Osdw = 
j2^f:2^kfx^l^^^^_^ Ogg^ ^ '£ss'^ss'i^-,si^+,s', and 
Otsc — '/'-,T''/^+,T (^"^ chose the 3^ = 1 component for 
convenience, the other components must be degenerate 
by su{2) spin symmetry). Calculation of the correlation 
functions associated with these operators is more subtle 
than in the spinless case due to the presence of the non- 
linear term in (jlOp . It was carried out in Ref. ^18,] using 
RG analysis. 

The resulting phase diagram, in the sense of diverging 
susceptibilities, is of the following structure: For positive 
backscattering gi± > fiows to 0+, 1, and Kp does 

not flow. To derive the power law decay of order param- 
eter correlations one can assume the fixed-point values 
of the Luttinger parameters. The result is 1/x^''^^ for 
CDW/SDW and l/a;^^"'+^ for SSC/TSC. This implies 
equally diverging susceptibilities of CDW and SDW for 
Kp < 1 and of SSC and TSC for Kp > 1. The degener- 
acy between order parameters is removed by logarithmic 
corrections favoring SDW and TSC [3]. In the opposite 
regime gi± < 0, the non linear term is relevant, opening 
a spin gap in the spectrum. Consequently SDW and TSC 
correlations are short ranged. The scaling exponents of 
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FIG. 3: Noise correlations of a spin-1/2 fermion Luttinger liq- 
uid with gi± < 0. Different rows correspond to different spin 
dependent correlations: (a,b,c) Gn, which in the spin gapped 
phase is sensitive only to CDW correlations; (d,e,f) Gn, sensi- 
tive only to the SSC correlations; (g,h,i) G = X^^^ Ga^ sensi- 
tive to both CDW and SSC. Columns correspond to different 
values of the Luttinger parameter Kp {Ka is set to zero in 
the spin gapped phase): (a, d, g) Kg, — 0.8; (b, e, h) Kp = 1; 
(c, f, i) Kp — 1.25. The noise correlations in G clearly show 
coexistence of the two (quasi) orders for 1/2 < A' < 2 (g,h,i). 

the remaining CDW and SSC correlations may be found 
by formally setting K^- ^ Q, giving R{x) ~ |a;|~'^*^ 
for CDW/SSC respectively. Hence there is a sizeable 
region ^ < Kp < 2, where both susceptibilities are diver- 
gent. 

The expressions for the spin-resolved noise correlation 
functions G||(q, g') and Gii{q,q') are identical to Eq. 
(I6|), with the exponents in Eq. ^ given by: 

g = (Kp + K;')/8 + {K, + K-')/8-l/2 
/iTT = {Kp-K;')/8 + {K,-K-^)/8 
hn = {Kp-K;')/8-{K,-K-^)/8 (11) 

As before, we note that this complicated function con- 
tains, along with other contributions, the equal time cor- 
relations of the four order parameters discussed above. 
For example G-f^ contains the TSC and a component of 
the CDW correlations, whereas contains the SDW, 
as well as components of both the singlet and triplet SC. 
Thus, G = X^ss' ^ss'i contains all the order parameter 
correlations equally. As in the spinless case we expect 
that when these correlations decay slowly enough they 
give the most singular contribution to Ggs' ■ Where they 
are not singular, other contributions may take the lead. 

Let us first consider the case gij. > taking the val- 
ues of the parameters at the fixed point: K^ 1 and 
9i± 0. We note, that in contrast with the spinless 
case, the equal time order parameter correlations decay 
too fast to give a singular contribution to Gss' ■ Exactly 
as in the spinless case we can obtain analytic expressions 
for the leading contributions to Gss' for Kp w 1 . Expand- 
ing to quadratic order in hss' we again find Gss'i<l,Q') 



given by Eq. ([8]) with hss' replacing h. This general 
behavior agrees with the result of numerical integration 
of Gss'{q,q'), although there may be corrections to the 
power of the singularities away from the line Kp = 1. 

In the opposite regime of negative backscattering one 
should evaluate the integral ([6]) with the exponents given 
by (fTTjl in the limit K^ — > 0. Taking, for example, G|f we 
see that hi-\ — > — oo effectively enforcing xi4 ± Xi2 = 0. 
The remaining integral over x = {xi2 — X3i)/2 then de- 
pends only on the combination 4(17 -|- /i^-f) + 2 = Kp. 

We find Gn sgn(l - Kp)\q - q'l^"-^ and G|| ~ 

sgn(l — 1/Kp)\q — g'l^/^"^^. These are precisely the ex- 
pected contributions from equal time CDW and SSC or- 
der parameter correlations. In the region 1/2 < Kp < 2 
where the susceptibilities of both order parameters di- 
verge, we find a divergence in the noise correlations due 
to only one of the orders and a sharp cusp due to the 
other incipient order. Specifically, for 1/2 < Kp < 1, we 
find a divergence in G||, and a cusp in G||, whereas for 
1 < Kp < 2, we find a cusp in G||, and a divergence in 

We also carry out a numerical integration of The 
results are displayed in Fig. ([3|). They are in perfect 
agreement with the analytic forms discussed above. Pres- 
ence of a spin gap is immediately seen by inspecting 
the different components of Gss' ■ Specifically, particle- 
particle correlations arise only in the ti channels, there- 
fore they correspond to SSG. Similarly particle-hole cor- 
relations are seen in the tt and not in the ti channel, 
consistent with CDW but not SDW quasi long range or- 
der. 

In summary, we propose that analysis of the spatial 
noise correlations in time of fiight images, may serve as 
a detailed probe of correlations in one dimensional sys- 
tems of ultra-cold fermionsQ. In contrast with Bragg 
scattering, which can only detect correlations in the 
particle-hole channel (CDW/SDW), such a probe would 
be equally sensitive to superfluid correlations in the 
particle-particle channel. The power-law singularities in 
the noise correlations are directly related to the singu- 
larities of the equal time order parameter correlations, 
hence they allow to extract Luttinger parameters from 
experiments. In addition we find new singularities in the 
noise correlation function, pointing to the ground state 
structure, even where the pure equal time order param- 
eter correlations are non singular. 
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